Effective slip over superhydrophobic surfaces in thin channels 



00 
O 

o 

< 



>> 
i 

q 

O 

CZ5 
>Y 

^ : 

Ph. 



> 

m 
m 

00 

o 

oo 
o 




Francois Feuillebois, 1 Martin Z. Bazant, 1,2 ' 3 and Olga I. Vinogradova 1 ' 4 ' 5 

} CiViJ5 [/Mi? 7636 and 7083, ESPCI, 10 rue Vauquelin, 75005 Paris, France 
2 Department of Mathematics, Massachusetts Institute of Technology, Cambridge, MA 02139 USA 
J Department of Mechanical Engineering, Stanford University, Stanford, CA 94305 USA 
* A.N. Frumkin Institute of Physical Chemistry and Electrochemistry, 
Russian Academy of Sciences, 31 Leninsky Prospect, 119991 Moscow, Russia 
S ITMC and DWI, RWTH Aachen, Pauwelsstr. 8, 52056 Aachen, Germany 
(Dated: August 10, 2008) 

Superhydrophobic surfaces reduce drag by combining hydrophobicity and roughness to trap gas 
bubbles in a micro- and nanoscopic texture. Recent work has focused on specific cases, such as 
striped grooves or arrays of pillars, with limited theoretical guidance. Here, we consider the ex- 
perimentally relevant limit of thin channels and obtain rigorous bounds on the effective slip length 
for any two-component (e.g. low-slip and high-slip) texture with given area fractions. Among all 
anisotropic textures, parallel stripes attain the largest (or smallest) possible slip in a straight, thin 
channel for parallel (or perpendicular) orientation with respect to the mean flow. For isotropic (e.g. 
chessboard or random) textures, the Hashin-Strikman conditions further constrain the effective slip. 
These results provide a framework for the rational design of superhydrophobic surfaces. 

PACS numbers: 83.50.Rp, 47.61.-k, 68.08.-p 



Introduction.- The design and fabrication of micro- 
and nanotextured surfaces have received much attention 
in recent years. It has also been recognized that a mod- 
ified surface profile can induce novel wetting properties 
of a solid, which could not be achieved without rough- 
ness [![. Depending on intcrfacial characteristics, the 
Wenzel state, where the liquid impregnates the surface, 
can enhance wettability, or the Cassie state, where the 
texture is filled with gas, can dramatically amplify hy- 
drophobicity 0. The remarkable mobility of liquids 
on such superhydrophobic surfaces renders them "self- 
cleaning" and causes droplets to roll (rather than slide) 
under gravity and rebound (rather than spread) upon 
impact. Beyond their fundamental interest, superhy- 
drophobic surfaces may revolutionize microfluidics 0, [j] , 
by reducing viscous drag in very thin channels and am- 
plifying transport phenomena [Bf and transverse flows Q . 

Reduced wall friction is associated with the break- 
down of the no-slip hypothesis. It has recently become 
clear that liquid slippage occurs at smooth hydropho- 
bic surfaces, as described by the Navier boundary con- 
dition 

@, H II v s = bdv/dz, where v s is the slip (tan- 
gential) velocity at the wall and the axis z is normal to 
the surface. A mechanism for dramatic friction reduc- 
tion involves a lubricating gas layer of thickness S with 
viscosity \x g much smaller than that of the liquid \x [Io| . 
so that 6 ~ 5(jLt//Lt ff — 1) 505 [ll[. This scenario al- 
lows to achieve slip length of only of a few tens of nm 
in case of smooth hydrophobic surfaces The pres- 

ence of a rough texture however stabilizes the gas layer, 
and by increasing its height S, the slip length may reach 
tens of /im over the gas regions. The composite nature 
of the texture, however, requires regions of lower slip (or 
no slip) in direct contact with the liquid, so the effective 
slip length of the surface b* (defined below) is reduced. 



FIG. 1: Sketch of a thin channel, where the gap width h is 
small compared with the texture characteristic length L. 



For anisotropic textures b* depends on the flow direction 
and is generally a tensor [l3l ]. Indeed, experimental stud- 
ies of flow past superhydrophobic surfaces suggest that 
b* does not exceed several /jm 1A\ and varies with the 
orientation of the wall texture relative to flow 11511. 



The quantitative understanding of liquid slippage past 
superhydrophobic surfaces is still challenging, and little 
theoretical guidance is available for the design of opti- 
mal textures. Some exact solutions are known for a flow 
on alternating (par allel or transverse) no-slip and per- 



il, 1?], 18 1 or transverse inhomogeneous 



feet slip stripes 
slip sectors [19( ■ Simplified scaling ex pres sions have been 
proposed for a geometry of pillars fo L |20| | , an d numerical 
approaches have also been followed [2ll , 22 , 23| . Never- 
theless, general principles to maximize or minimize the 
effective slip have not yet been established, even in the 
simple (but experimentally relevant) lubrication limit, 
where the implication of slip is the most pronounced [Io| . 

In this Letter, we propose a systematic approach to 
optimize the effective slip length of a superhydropho- 
bic surface in a thin channel, based on the theory of 
heterogeneous porous materials 2J, |25j . We derive rig- 
orous bounds on the effective slip length for arbitrary 
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anisotropic or isotropic textures, depending only on the 
area fractions and local (any) slip lengths of the high-slip 
and low-slip regions. In some cases, the bounds are close 
enough to render detailed calculations unnecessary, and 
in others the theory provides optimal textures which at- 
tain the bounds (notably the maximum possible effective 
slip). Our theory also predicts b* in certain geometries 
without requiring any calculation. 

Model and analysis.- We consider pressure-driven 
flow of a viscous fluid between two textured parallel 
plates ("+" and "-") separated by h, as sketched in FigfT] 
Motivated by superhydrophobic surfaces in the Cassie 
state, we assume flat interfaces (such as idealization has 
been used in most of previous studies 3, 21 1 and corre- 
sponds to a minimum dissipation in the system 18,[26j]) 
characterized by spatially varying slip lengths b + (x,y) 
and b~(x,y). Our analysis is based on the lubrication 
(or Hcle-Shaw) limit of a thin channel, where the texture 
varies over a scale L 3> h; the flow profile is then locally 
parabolic at any position. 

To evaluate the effective slip length, we calculate the 
velocity profile and integrate it across the channel to ob- 
tain the depth-averaged velocity U in terms of the pres- 
sure gradient Vp along the plates. As usual for the Hclc- 
Shaw cell, the result may be written as a Darcy law 



U 



k(x,y) 



/' 



Vp, 



(1) 



where we obtain the permeability 



3Q3+ + /T +4/3+/3~) 
1 + (3+ + 



in terms of the normalized slip lengths /3 + = b + (x,y)/h 
and (3~ = b~(x,y)/h. The permeability is maximized 
with two equal surfaces, (3 + = (3~ = f3(x,y), so we 
consider this case (II) with the goal of minimizing drag. 
We also consider the case (I) of one no-slip wall ((3 + = 
f3(x,y); f3~ = 0), which is relevant for various setups, 
where the alignment of opposite textures is inconvenient 
or difficult. The permeability then takes the form: 



3(3(x,y)/[l 
6/3{x,y) 



f3(x,y)} 



case (I) 
case (II) 



(2) 

In general, the slip length may also vary locally with 
orientation, so that b(x, y) becomes a second-rank tensor 
h(x,y), from which a tensorial permeability h(x,y) can 



be derived [13] 



The slip length b(x,y) (or h(x,y)) varies on the mi- 
croscale L ^> h, but we are interested in properties of the 
flow at the macroscale. A natural definition of the effec- 
tive slip length is based on a hypothetical uniform chan- 
nel with the same effective permeability. First, we aver- 
age (P) over the texture (denoted by (.)) at a mesoscale 



that is smaller than the macroscale, but much larger than 
L, to obtain 

(U) =-hk(x,y)Vp) = — -(V P ) 
fj, /i 

where in the last step we introduce the effective perme- 
ability k*, which is generally a tensor, even if k(x,y) is 
locally isotropic. Only with an isotropic structure at the 
mesoscale docs it become a scalar k* . This definition 
is subject to the boundary condition of a uniform pres- 
sure gradient VP applied at the macroscale, which must 
equal the average pressu re g radient, (Vp) = VP, since 
the pressure is harmonic |24[. 

By analogy with |J2J), we define the effective slip length 
in terms of the effective permeability: 




case (I) 
case (II) 



(3) 



where the principal (eigen)directions j = 1,2 of k* corre- 
spond with those of 3* = h*/h, where b* is the effective 
slip length tensor [131 ] . 

Motivated again by superhydrophobic surfaces in the 
Cassie state, we assume 6(x, y) switches between two val- 
ues, bi and 62, associated with permeabilities fci, &2 from 
([2]) , for regions (or "phases" ) of liquid-solid and liquid-gas 
interfaces, respectively. Let <pi and 4>2 be the area frac- 
tions of the two phases with <j}\ + (f>2 = 1. We make no 
further assumptions in deriving bounds on the effective 
slip length (3* in a principal direction (without transverse 
flow), aside from distinguishing between anisotopic and 
isotropic textures. 

Anisotropic textures. In the general case of an 
orientation-dependent texture (k* 7^ k*T), the Wiener 
bounds apply for the effective permeability in a given di- 
rection [25(: k 1 - < k* < fell. The lower bound k 1 - can be 
attained by parallel stripes perpendicular to the pressure 
gradient: k 1 - = (0i/fci + 4>2/k2)~ 1 - The bound fc" can 
also be attained, by stripes parallel to the pressure gradi- 
ent: fell = 4>\ki +4>2k2- Physically, these special textures 
act like resistors in series and in parallel, respectively. 

Using ([2|) and ([3]), the corresponding bounds for the 
effective slip length are 

= — < a < = — case (I) (4a) 

l + 4(/3) -' ~ l + {(3) 

((3) + 60^2 



1 + 6(f3) 



< (3* < ((3) case (II) (4b) 



where 



(13) = cMi + 02/?2 and 0) = cfofc + fafa (5) 

are the average slip length and average transposed slip 
length, respectively. Using parameters for typical super- 
hydrophobic surfaces, these bounds are plotted versus 
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FIG. 2: (a) Bounds on the (normalized) superhydrophobic 
slip length fa" I fa versus the liquid-gas area fraction </>2, as- 
suming no slip fa = and high-slip fa = 5 on the liquid-solid 
and liquid-gas interfaces, respectively. Bottom: zoom of top 
figure. Dashed and solid lines correspond to cases (I) and (II) 
or one or two superhydrophobic surfaces, respectively. In each 
case, curves from top to bottom represent: the upper bound 
for anisotropic, upper bound for isotropic, lower bound for 
isotropic, lower bound for anisotropic textures. The value of 
fa for the chessboard or the isotropic Schulgasser structure 
sketched in Fig. [3] is also shown (square: case (I), circle: case 
(II)). (b) The same bounds plotted versus the slip length fa 
for 4> 2 = 0.9. 



the liquid-gas area fraction fa in Fig. [2ja) and versus 
the liquid-gas slip length fa in Fig. EJb) . In case (I) the 
bounds are fairly close (especially when fa is large), so 
the theory provides a good sense of the possible effec- 
tive slip of any texture, based only on the area fractions 
and local slip lengths. In case (II) the difference between 
the upper and lower bounds is larger and grows quickly 
with fa. In either case, however, the texture attaining 
the upper (lower) bound corresponds to stripes oriented 
parallel (transverse) to the pressure gradient [27 1. 



Isotropic textures.- Consider now any isotropic 
structure, without a preferred direction (k* = k*T). If 
the only knowledge about the two-phase texture is fa, fa, 
then the Hashin-Shtrikman (HS) bounds apply for the ef- 
fective permeability, k^ s < k* < k^ s , where (assuming 
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FIG. 3: Special textures arising in the theory: (a) stripes, 
which attain the Wiener bounds of maximal and minimal ef- 
fective slip, if oriented parallel or perpendicular to the applied 
pressure gradient, respectively; (b) the Hashin-Shtrikman 
fractal pattern of circles, which attains the maximal slip 
among all isotropic textures; and (c) the Schulgasser texture, 
whose effective slip follows from the phase-interchange theo- 
rem. 



fa < fa without loss of generality) 

1-12 

'^HS 



k L 
K HS 



(k) 



(k) + fa 



[kf 



(k) + fa 



with [k] = fa — ki and using the same notation as in 
J5|). Using (|2|) and ([3]), the corresponding bounds for the 
effective slip length are obtained in a form similar to ([4]): 

(P) + f {fa) fa fa <(] * < (fa" + f (fa) fa fa (6) 



i + /GW> 



1 + /GW) 



where 



/(/?) 



(5 + 3/3)/(2- 
3/(1 + 3/3) 



5 fa 



case (I) 
case (II) 



(7) 



The HS bounds © are plotted in Fig. [2] in the same way 
as the Wiener bounds (g]) and behave similarly, aside 
from being closer and confined between them. However, 
it turns out that isotropy does not dramatically reduce 
(enhance) the maximum (minimum) effective slip in a 
thin channel, especially in the configuration with two su- 
perhydrophobic surfaces, case (II). 

The upper bound for isotropic textures can be attained 
by a fractal pattern of nested circular patches [25| as 
shown in Fig [3] It is interesting to note that similar 
patterns might be expected for a random (and some- 
times fractal) nanobubble coating [28|. However, it is 
not necessary to deal with fractal surfaces: some periodic 
honeycomb- like structures also attain the bound [i^] . 

Finally, we use phase interchange results 25[ to obtain 
the effective slip length without any calculations, for a 
special class of isotropic textures. For a medium that 
is invariant by a 7r/2 rotation followed by a phase inter- 
change, a classical result follows: k* = \ / k 1 k 2 - Examples 
of such media are the chessboard and the Schulgasser 
texture, sketched in Fig. [3]Jc). The effective (3* is then 
easily obtained from ([3]), although the corresponding val- 
ues, shown in Fig. (2[ are far from the HS upper bound. 

Concluding remarks: design strategies- We close 
by proposing some guidelines for the design of thin su- 
perhydrophobic microchannels, which maximize effective 
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slippage, e.g. for lab-on-a-chip applications. We as- 
sume a principal direction of the texture is aligned with 
the side walls, since this is typically the fastest orien- 
tation. (Tilted textures also complicate analysis, since 
the constraint of no transverse flow, (U) y = 0, induces a 
transverse pressure gradient, (Vp) y = —(k* x /k*)(Vp} x , 
which in turn affects the mean forward flow, (U) x — 
det(\t*)/k* y 0.) For simplicity, we also restrict 
now to the case j3\ = of no-slip support structures. 
It has been predicted for thick (L <C h) cylindrical [l6| 
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and planar channels [17|, l2l( that the longitudinal stripe 



configuration has larger effective slip than the transverse 
one. For a thin channel [L 3> h), we can now draw the 
more general conclusion that longitudinal (transverse) 
stripes provide the largest (smallest) possible slip that 
can be achieved by any texture. Interestingly, this in 
contrast to a prediction for thick channels, where an ar- 
ray of pillars in the limit <f>2 — ► 1 has larger slip than 
longitudinal stripes [2tj| . 

We have shown that the key parameter determining 
effective slip is the area fraction of solid, <f>i, in contact 
with the liquid. If this is very small (or <p2 — > 1), for all 
textures the effective slip tends to a maximum, /3* — > /52- 
In this limit, the microchannel produces a kind of su- 
perfluidity, with plug-like flow. However, even a very 
small (f>i is enough to reduce the effective slip signifi- 
cantly since in this limit (except an upper limit for case 
(II), where 2 — 0*/02 — 0) we have the asymptotic scal- 
ing — 0* 102 oc 02<fri- It is interesting that in case 
of perfect slip over the gas areas, 0* scales as oc <t>2/4>ii 
which is similar to an earlier result for a thick cylinder 
with transverse stripes For thin channels, we see 

now that this result is very general and is valid for any 
texture (and likely any channel geometry) with perfect 
slip patterns, representing "obstacles" to the flow. We 
thus conclude that in many situations, maximizing (32 
is not nearly as important as optimizing the texture to 
achieve large effective slip. 

Finally, we have demonstrated that for all slip lengths 
and all fractions the largest possible (3* is equal to the 
area-averaged slip length ((3), attained by longitudinal 
superhydrophobic stripes. However, for all textures the 
effective slip nearly coincide with the average, provided 
02 is small (or, more generally, 02 — f3i is small). Al- 
though this limit is less important for pressure-driven 
microfluidics, it may have relevance for amplifying trans- 
port phenomena 

In summary, we have connected the problem of effec- 
tive slip over superhydrophobic surfaces in thin channels 
with the classical subject of conduction in heterogeneous 
media. This has allowed us to obtain rigorous bounds on 
slip for arbitrary textures and to obtain the slip in some 
cases without any calculations. Our results can be used 
to guide the design of superhydrophobic surfaces for thin 
micro- or nano-channcls (where slip is most important), 
and some principles may hold for thick channels as well. 
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